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SAFETY STOCK DETERMINATION 

The invention relates to a method for assigning to a set of products a set of corresponding 
continuous densities of demand distributions. More particularly it relates to a method for 
determining safety stock levels for several products and for controlling the corresponding 
5 stock level. 

TECHNICAL FIELD AND BACKGROUND OF THE INVENTION 

Warehouses are used as stock supply of on-sale items because of uncertainties of demand or 
because some suppliers only deliver large lot sizes. Apart from figuring out what product 
has to be ordered when and at what quantity, one of the most important problems in 

10 warehousing is to set the size of the safety stock correctly. The safety stock is supposed to 
cover the uncertainty of demand during the inbound lead time. If the safety stock level is too 
high then more financial resources are bound to the inventory than necessary. If, however, 
the safety stock level is too low, the supplier will more often run into out-of-stock situations 
and thus will not be able to satisfy the demand knd thereby fail to meet the desired 

15 serviceability. 

Most of the software tools used in supply chain management and warehousing determine 
the level of safety stock for each of the products under the assumption that the daily 
demands are normally distributed or by simply setting a fixed safety stock level- However, 
the problem with this assumption is that in many jcases the demand is far from being 
20 normally distributed. In particular, this is true for certain kinds of warehouses like for spare 
parts. In the end this results in incorrectly calculated safety stock levels, i.e. either the 
achieved service levels end up being too high and therefore the inventory level could have 
been reduced or the achieved service levels end up being too low, thereby failing to meet 
the desired serviceability. In both cases additional costs are generated. 

25 From Kohonen, Teuvo: "Self-Organizing Maps", 3rd -edition, Springer Ser. in Ihfomiation 
Sciences 30, Berlin 2001, a clustering algorithm is known. Additional information thereto 
can be found in Kohonen, Teuvo, 'The self-organizing map'\ Neurocomputing 21, No. 1-3, 
1-6 (1998). 

I 

i 
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In KuJiback, Leibler: r, On Information and Sufficiency", in Annals of Mathematical 
Statistics 22 (1951), pp 79-76, the Kullback-Leibler distance is described. 

In Scarf, "The Optimality of (s,S) Policies in the Dynamic Inventory Problem", in Arrow, 
Suppes (Ed.), Mathematical Methods in the Social Science", Stanford University Press, 
5 Stanford I960" the (s,S) replenishment method is described. More about that can be found 
in Zheng, Y.-S. "A simple Proof for the Optimality of (s,S) Policies in Infinite Horizon 
Inventory Systems", Journal of Applied Probability 28, 802-810 (1991) 

In Schneider, H. t "Effect of Service-Levels on Order-Points or Order-Levels in Inventory 
Models", International Journal of Productions Research 19, 615-631 (1981) service level 
10 definitions are described. 

OBJECT AND ADVANTAGES OF THE INVENTION 

According to an aspect of the invention as set forth in the claims a method for assigning to a 
set of products a set of corresponding continuous densities of demand distribution is 
15 provided. 

The proposed method has the advantage of offering a safety stock level with higher 
precision, i.e. offering a product stock controller a means to get closer to his target to 
minimize the overall safety stock value under the condition of maintaining a given service 
level. The precision is particularly achieved utilizing a clustering step that on one hand 

20 reduces the number of determinations in a subsequent selection step and that on the other 
hand is able to introduce a cross-product statistical data enhancement effect, i.e. improve 
the demand statistics especially for sporadically demanded products. This allows to use 
more density models in the selection step and thereby achieve better precision in that step. 
Using the same number of density models for direct comparison with the products* densities 

25 would lead to an excessive uneconomic determination workload. 

According to another aspect of the invention a method is provided that further encompasses 
a safety stock determination wherein with a given service level from the continuous demand 
densities a product safety stock level is determined. 

According to yet another aspect of the invention the method can be enhanced by a product 
30 order step wherein via an output device an order for a product is initiated towards a 
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product-providing entity, e.g. via a communication medium like a communication network, 
in the event that a stock level of that product falls below a reorder point. Thereby an 
automated product stock control is achievable. 

The method can further comprise a clustering correction step wherein for each product a 
5 distance between its continuous density and its normalized discrete demand density, is 
determined and wherein for those of said products whose distance exceeds a predetermined 
threshold value, in a substitute-selection step for the corresponding normalized discrete 
demand density out of said predetermined set of continuous model densities a substitute 
continuous density is selected, which is defined for that product as its continuous density. 
10 That cluster correction step is advantageous because it is able to perform a distance check 
between the determined continuous densities and the normalized discrete demand densities 
and perform a correction should the distance be found too large. This improves the 
precision of the resulting assigned continuous densities and thereby reduces the error in the 
final safety stock levels. 

15 In the conversion step the demand time series can be substantially segregated from demand 
trend and periodic demand. Segregating the demand time series from demand trend and 
periodic demand is another way of increasing precision, since the parrs that are removed 
would lead to eirors in the resulting assigned continuous densities and hence in the final 
safety stock levels. 

20 In the normalization step normalization information can be maintained for the adjustment 
step. Maintaining the normalization information for the adjustment step facilitates the 
method since that information can be directly used there and need not be constructed in a 
different manner. 

If for each of the clusters an initial cluster representative is provided, and if Lhe clustering 
25 step comprises for each of the normalized discrete demand densities a distance 
determination step wherein for each said initial cluster representative its distance from the 
normalized discrete demand density is determined, and a modification step wherein the one 
initial cluster representative is modified towards the normalized discrete demand density to 
which it has the smallest distance, said modified cluster representative being thereafter 
30 defined as new initial cluster representative. Using a difference estimation step and a 
modification step in the clustering step provides advantageous because this method is a 
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quick and precise way of clustering. Furthermore the Kohonen algorithm can be utilized 
therefor. 

If the distance determination step and the modification step are repeated one or more times 
for each of the normalized discrete demand densities, starting with the initial cluster 
5 representatives that are present after the previous last modification step, this provides the 
advantage that the clustering result will be preciser which again leads to less errors in the 
resulting assigned continuous densities and hence in the final safety stock levels. 

Selecting the initial cluster representatives that are present after the last modification step in 
the clustering step as the cluster-representative discrete density brings the advantage that 
10 since the initial cluster representatives are already present, no additional steps for the 
selection have to be performed. This hence reduces the rime the method takes. 

If in the clustering step a said cluster representative is substantially not modified when its 
distance to those normalized discrete demand densities that belong to the corresponding 
cluster is below a predetermined value, unnecessary representative movements are avoided, 
15 which on one hand reduces the clustering time and on the other hand avoids an exacerbation 
of the resulting clusters. Adding a condition that this cluster need have a number of 
thereto-belonging demand densities, also called members, that is bigger than a 
predetermined number, would introduce the advantage that only clusters of a predetermined 
minimum size are handled this way, again improving the clustering result. 

20 If the predetermined set of continuous model densities contains the model densities with its 
parameters in an undetermined form and the parameters for the clusters, are determined in 
the selection step, this brings in the advantage that the selection step becomes less 
storage-intensive since less model densities need be stored. 

If in the normalization step the demand value for 0 pieces of the products in the discrete 
25 demand densities is reseated inversely to the demand values for more than 0 pieces of the 
same product, a preciser clustering for products with sporadic demand, also referred to as 
slow-moving items, is the advantageous result. 

If in the selection step the demand value for 0 pieces of the products is suppressed, the 
selection result is less distorted by that demand value, which improves the resulting 
30 selection. 
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If the method further comprises a safety stock determination step wherein with a given 
service level from the continuous demand density for one or more of the products a product 
safety stock level is determined, a stock controller gets preciser information about how 
many pieces of the product should be kept on stock in order to guarantee the service leveL 

5 Performing a product order step wherein via an output device an order for a said product is 
initiated towards a product-providing entity, in the event that a current stock level, of that 
product falls below a minimum quantity that is the sum of the average of the discrete 
demand density for said product during the product's lead time, and the determined safety 
stock level for said product provides advantageous since it realizes an automated stock 
10 control for the products, thereby decreasing the probability of out-of-stock situations. 

The above methods are implementable in form of program code means for performing the 
method. 

SUMMARY OF THE INVENTION 

The invention is directed to a method for assigning to a set of products a set of 
15 corresponding continuous demand densities. In a conversion step for each product its 
demand time series is converted into a discrete demand density. In a normalization step the 
discrete demand densities are transformed into normalized discrete demand densities. In a 
clustering step each of the normalized discrete demand densities is assigned to a cluster and 
for each cluster a cluster-representative discrete density is determined. In a selection step for 
20 each cluster-representative discrete density out of a predetermined set of continuous model 
densities a cluster-representative continuous density is selected. In a 
parameter-determination step for each product for its cluster-representative continuous 
density product-individual density parameters are determined under use of which for each 
product a continuous density is determined. In an adjustment step for each product the 
25 continuous density is adjusted to form therefrom Che continuous demand density. 

DESCRIPTION OF Tip: DRAWINGS 

Examples of the invention are depicted in the drawings and described in detail below by 
way of example. It is shown in 
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Fig. I a block diagram, specifying the different steps and results of a density assignment 
process. 

Fig. 2 a schematic view of a computer in which from given product-specific input 
parameters a stock-alert signal or an automatic product order are issued. 

All the figures are for sake of clarity not shown in real dimensions, nor are the relations 
between the dimensions shown in a realistic scale. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

In the following, the various exemplary embodiments of the invention are described. 

For sake of simplicity, in the following, when referring to a "density 1 ', the density of the 
demand distribution is meant. When the index "i" is used herein for a specific reference 
indicator this is meant as a variable that could take any value of the specific values used 
elsewhere in this description for the same reference indicator. 

In figure 1 a block diagram, specifying the different steps and results of a density 
assignment process is shown. 

For each of a multitude of products Pi, P 2> Ps> P4, a corresponding demand time series H(P,), 
H(Pa)» H(p3>, H(P*) is given. The demand time series can be a historical demand data, e.g. 
the demand development for the respective product P,, P 2> P3, P4 watched over a specific 
period of time. The time is herein quantified into time units, e.g. days. The demand is herein 
quantified into demand units, i.e. the smallest measured demand unit, to allow comparison 
of different product types. E,g. for screws a package of 1000 pieces can be regarded as 
smallest traded and demanded unit, whereas for orange juice, 1 liter will be regarded as 
smallest traded and demanded unit. These demand units will in the following be referred to 
as demand quantities. 

In a conversion step 10 for each product P u P2, Pj, P4 its demand time series H(Pj), H(P 2 ), 
H(P 3 ) ? H(p4> is converted into a corresponding discrete demand density D PJ a, Dkm, D p ^, D P4 d 
by a binning process that defines a finite set of demand bins into which the actual demands 
are put. Such a bin could e.g. provide for a demand from 1 to 5 pieces belonging into a first 
bin, 6 to 10 pieces belonging to a second bin* a.s.o. Each of these discrete demand densities 
Dpi* D^ 2d , Dp3d, Dp4d hence shows for its corresponding product P t , P 2 , Pj, P 4 . how often a 
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specific amount of that product P. u P 2 > P 3t P 4 is demanded. This discrete demand density 
Dpu, Dpm, Dp*>, D P4 fl is free from demand trend and periodic demand. The purging can be 
done in the conversion step 10, but can also be done before, i.e. the demand time series 
H(P,), H(P 2 ), H(P 3 ), H(P 4 ) can also be already provided in a form that is free from demand 
5 trends and periodicities, "Free" herein means substantially free, since the purging in reality 
can not guarantee a 100% removal of trend and periodicity. Thus in the discrete demand 
density D intS , Dpw, Dpw> Dp*! substantially the stochastical component remains. Since 
demand trend and periodic demand are those demand components that can be more easily 
predicted, the safety stock need not cover it. The remaining, stochastical component is 
10 hence the one for which the safety stock is provided for the corresponding product P tt P 2 , 

As next step follows a normalization step 20 wherein the discrete demand densities Dpi d , 
Dp2d, Dp3d, Dp4d are transformed into normalized discrete demand densities D^Vid. D N P 2d» 
D*W D N i*d. The normalization provides for a scaling to a norm value for the average and a 

15 norm value for the variance. The typical norm value for the average is 0 and the typical 
norm value for the variance is 1. The normalization to the norm value for the average 
effects a shifting by a shifting-factor s*, whereas the normalization to the riorm value for the 
variance effects a scaling by a scaling-factor Sy- These two factors So, Sv can be maintained 
for each normalized discrete demand density D N P id, D N p2d, D^W respectively the 

20 corresponding product Pi, P 2 , P 3 , P* as adjustment factor set v tt] , v„2, v s & t v** for a later 
adjustment step 70. The normalization has the effect that in the subsequent step, the 
densities D N ?id, E^ra, D N P3* D^W that are compared have a common scale which effects a 
preciser relative positioning of the densities D^a, D N p 2 d» D"™, D^W with respect to each 
other. More precisely, the distances idcCX, D*™), as explained below, will be more accurate 
* 25 which leads to a fairer clustering, iie. those densities D N P i d , D*W D**^, D"^ which have 
the same shape will eventually be assigned to the same cluster. 

As next step follows a clustering step 30 wherein each of the normalized discrete demand 
densities D N Hd, D'W D^W D K pa<t is assigned to a cluster Ci, C 2 and wherein for each cluster 
C. C 2 a cluster-representative is determined which is a cluster-representative discrete 
30 density Ric. The number of clusters Ci. C 2 is predefined and can be determinable for this 
method. 
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The clustering step 30 comprises for each of the inonnalized discrete demand densities 
D N Pi d , D 1 *™, D N P3d, D N wd two steps, a distance determination step and a modification step. 
First for each of the clusters C Jp Cz an initial cluster representative Ri Cl also referred to as 
initial cluster representative Xu X2, is provided. The initial cluster representative Xu X2 
5 need not be a normalized density since in the subsequent steps the quality of being 
normalized does not play a role. Normalization is here only of value for the relative 
positions of the discrete demand densities D N nd, D^W D N PW > D^W 

i 

: i 

In the distance determination step, for all initial cluster representatives Xi, X 2 their distance 
dcCXj, D N Rd) from one of the normalized discrete demand densities D N pi d » DVw* D%, D"^ 
10 is determined. The one normalized; discrete demand density D N P t d , l^W D^W D N P4d can be 
selected according to any selection mechanism, bei it a simple sequential scheme, or a 
random selection. 

As distance herein the KuUback-Leibler distance c^n be used, by defining each cluster 
representative X>, X2 and each normalized discrete demand density I^W D N ra , D N P3d, D n p*i 

15 as a point in an n-dimensional Kullback-Leibler space, wherein n is the number of quantity 
bins used in the conversion step' 10. The position! of the normalized discrete demand 
densities D N P!4 , D N P2d, D^a, D N P4 <i in the Kullback-Leibler space depends on their average 
and variance, which due to the normalization step 20|has been unified, thereby leading to a 
common determination basis and acorrect relative positioning within the Kullback-Leibler 

20 space. The positive effect of the normalization step 20 is present for any type of clustering 
involving distance-determination, j. ! 

The start position of the initial cluster representatives Xi, X* can be preselected, but it is 
found that performing a number of iterations, as is described below, the start position 

! i 

becomes practically irrelevant. Therfore a start position of the initial cluster representatives 
25 Xi, X 2 can be selected, that is independent from the tyjpe, number or density of the products 
Pi, P2, P3, Pd. The method can beiidesigned to renounce a possibility to determine start 

positions and thereby require less user interaction by ajstock controller. 

j 

Next follows a modification step. this step among I the found distances dc(Xi, D N W d) the 
smallest distance and the corresponding pair of initial cluster representative X t and 
30 normalized discrete demand density D N * a is determined. Then that initial cluster 
representative Xi is moved, i.e. modified, towards thatj normalized discrete demand density 



EMPFANGSZEIT 29. MAR, 17:00 AUSORUCKSZE I T 29. MAR. 17:14 



FR 17:00 FAX +41 1 724 S9 51 



IBM ZURICH IPD 



10 



IS 



20 



25 



30 



-9- 



EP01 PATENTS 



CH 9-2001-0064 



@016 



D N rw> if the distance dc(Xi, D N wd) Is not already below a predetermined maximum distsunce 
level, e.g. 0. The Kohonen clustering algorithm pan be u$ed to perform that function, 
wherein the cluster representatives Xi will be what it the Kohonen algorithm is referred to 
as centers. 



The distance determination step and the modification step are performed for each of the 
normalized discrete demand densities D N nd, Wru* D^wa, D*W In the case of only a single 
run of the distance determination step and the modification step, it is useful to guarantee 
that every normalized discrete demand density D N i>td,;D K ?2d» D*™, DVa is once handled this 
way In case of multiple iterations, a random selectidn is deemed to satisfy this criterion by 
effect of its statistical uniform distribution. In that case, a tracking of the use of each 
normalized discrete demand density D*W D K P2d, E>*W D N p« is not necessary. 

The clustering step 30 hence has the effect that the clusters d, C 2 and their representatives 
Xt, X 2 are defined and that furthermore it is definled which of the normalized discrete 
demand densities JD 1 *™, D N P2d, D^Vm, D M i*d belongs to Which of the clusters Ci. C 2 . 

The clustering step 30 can be repeated in that in the distance determination step for each 
modified cluster representative Xi, X 2 its distance dc(Xj, E^rw) from the normalized discrete 
demand densities D N Pld , D",^, E^rea* E>"p*i is deteniiiied, and in the modification step each 

, i 

modified clusrer representative X ly & 2 is modified towards the normalized discrete demand 
density D*™, E^pm, D^pm. D"™ jo which its distance dc(Xj, I^W) is smallest. This 
repetition, also called iteration, can be done several times, referred to as £poque in the 
Kohonen algorithm. Thereby the accuracy of the clustering step 30 is increased This means 
that the final clusters d ( C 2 are more precisely defined, meaning that the modified cluster 
representatives Xi> X 2 lie in a closer distance to the; normalized discrete demand density 
D N j>, d , D N P2d, D N P3d, D N «d that belong to its cluster G, C*. 

It has been found that irrespective of the number of pr Dducts Pi, P 2> P 3 . P* and even itregard 
of the density types thereof, a number of between 30 and 50 iterations are sufficient to 
achieve a sufficient stability. Generally it can hence be stated that it is advantageous to 
perform more than 30 iterations, siice with the number of iterations the execution time 
rises, it can also be stated that less then 50 iterations aj e advantageous. As stability criterion 
Ae average distance dc(Xi, between the modifi ?d cluster representatives X,, X 2 and 

the normalized discrete demand density T> H pn, D N iw W™> D N P4 d that belong to the 
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respective cluster Ci. C 2 , can be use d. Also an automktism can be used that stops performing 
further iterations if a certain stability of the sample locations has been reached, i.e. the 
average distance dc(Xi, D N Rd) has fallen below a predetermined maximum distance value: 

The distance determination step can be executed under use of different distance 
determinators, thereby providing several values for tit e distances d c (Xi, D N «d). In such a case 
the modification step is executed using the smallest ralue of the distance dc(Xj, D^W. The 
use of differenr distance determiii^tors allows to avoid situations in which the distance 



not correctly reflect the true relation 



d c (Xj, D N pj d ) that is calculated andj determined does 
between the modified cluster representative Xt, X 2 and the normalized discrete demand 
densities J^W D N ?2d, D>W D*W : The different distance determinators typically perform 
different distance-determination alg srithms that henc? also result in different values for the 
distances dc(Xi t D n pw). It is assumat le that the distanc e determinator with the smallest value 



of the distance d c (Xi> D N ™) is the 



nost accurate ana useful. The terra "distance" is to be 



understood as a metric in the mathematical sense. 



Finally a cluster-representative disci *te density Rtc is determined for each cluster Ci. C 2 . The 
simplest selection is to select the m xiified cluster representatives X u X 2 that exist after the 
last modification step as the cluster -representative d screte density R«:> This selection has 
the advantage that no additional algorithm for selecting the cluster-representative discrete 
density Ric has to be performed.; It is also assure iable that after running through the 
modification step, in particular aftei several iteration!, the modified cluster representatives 
Xu X 2 are the best possible choice fior the cluster-representative discrete density Rpc- Since 
the cluster-representative discrete density Re is only i ised for a calculation step, it need not 
be any of the nonnalized discrete de nand densities DM P3 d, r> N Md„ D^W The clustering 
step 30 has a very valuable advantage in that by the clustering products Pj, P 2 * P4 are 
grouped that not necessarily have an r correlation to each other. In particular this means that 
those products P,, P 2 , Ps, P* which h ive only few demand occurrences in the period of time 
used for the demand time series H(l ,), H(P 2 ), H(P 3 ), H(P4), in the following referred to as 
slow-moving items, profit from thos > products P,, P 2 . > 3 . P* tbat have a similar normalized 
discrete demand density D N ™, D N * a , D>W D N P4 d, tut which have a higher number of 
demand occurrences in the same j>eriod of time. (This increases the accuracy of the 
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subsequent selection step 40 fo) 
occurrences, 

I 

The purpose of clustering the dam is hence twofoliL First, it enhances the quality of the 
succeeding steps by enlarging the I imount of data that is being analyzed. This is especially 
true for slow-moving items becaus t those have only little statistical data and the clustering 
effects the use of statistical data 6 f all products Pi, P 2 , P 3 , P 4 in that cluster. Second, the 
subsequent steps need only be perfi rmed on the cluster's representatives X, X 2 . 

The clustering step 30 can be furtl er improved by j ntroducing an exceptional rule: In the 



is substantially not modified when its 
densities D*™, D*^ D n pw, D*** that 
listance value. By giving the distance 



clustering step 30, the cluster reprjj sentative X u X 2 
10 distance dc(Xi, Pi) to the normalized discrete demanij 
belong its cluster Q, C2, is below ik predetermined 

value for the distance dc(Xj> I?*m)* Jeing a maximum value, a cluster representative X,, X2 

ij 

is substantially prevented from beit g moved away fiom these normalized discrete demand 
densities D N Pld , D*W D N rad , D*™/' o which it is so < lose. This reflects the idea that in the 
15 case one of the cluster representati *es X u X2 is veiy close, Le. the distance dc(X, D N Kd ) 

1 , 

being below the predetermined di$ ance value, to me corresponding normalized discrete 
demand densities D*™, I> N m 7 D>] 9d , D 1 "**, then *ie cluster representative X, is to be 
regarded as a good representative ari i hence shall not be modified, i.e. not by some effect of 
the clustering step 30 be moved awj y from its good p ssition. This cluster representative X r , 
20 X 2 is practically frozen to its positj >n, avoiding unwwted exacerbation of its position and 
even reducing the determination tii ie in the modific rtion step. Additionally the condition 
that the cluster Q, C2 should live a minimuifci size, i.e. a minhrrwmv - number of 
thereto-belonging normalized discrae demand densities D*™, D^pm, D N i>3d, D N M d, may be 
added, to ensure that the cluster frei ring i$ only perfe rmed for clusters Ci, C 2 of significant 
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, P2, P3, P4 with only few demand 



size. 



30 continuous model densities D^mc <9 ie that fits besf 



As next step follows a selection si sp 40 wherein ftr each cluster-repr^entative discrete 
density Rtc out of a predetermii td set L D of < ontinuous model densities D*W a 
cluster-representative continuous dc tsity D m c c is sel&ted. This is hence the step that goes 
from the discrete densities to the | :ontinuous densS ics. The task is to find among the 

the respective cluster-representative 
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discrete density Rrc* The set L D should comprise mc c el densities of different basic types like 



Gauss, Gamma, Poisson, Weibull, 



- 12 



Logistic, and oth i r parametrized density families. 



The predetermined set Lp of corfcnuous model censities D m Mc here contains the model 
densities D'V with their parameter p in an undeterp lined form. This has the advantage that 
the set need not contain the model densities D% c in a lot of possible parametrically 
determined forms which would! unnecessarily increase the need for storage and 



4p 

i! 



determination time for going through all of them. Ir s tead, in the selection step 40, the mere 



shape of the model densities D m Mc is compared to tl t 



Rfc to identify a best fit. An algorijthm suitable for 
moments or Maximum Likelihood [Estimators (MLi 
it is possible to have some of the parameters p 
determined once they are selected for each cluster- 
is in particular useful for shape parameters. 
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cluster-representative discrete density 
this operation would be the method of 
. In the course of the selection step 40 
yf the model densities D^W already 
resentative discrete density Ric- This 



ispj 



After the selection step 40 follows a pararneter-det mnination step 50 wherein for each 
product Pi, P 2 , P 3 , P4 for [its cluster-repi4j entative continuous density D m Cc 



product-individual density parameters p(Pi) are del 2 rmined under use of which for each 



D^k, D"W is determined. This step 



product P u P2. P31 P4 a continuous <lensity D^pu, D m i 2t> 
hence leads back from the clusterw^e handling to th productwise handling 

The cluster-representative continuojis density D^cc tttit 
O. C 2 is now for each product Pi,jjP2. P3, P4 that hi 



compared to the corresponding normalized discrete 



D N P4d and by means of that comparison the product- j idividual density parameters p(Pi) for 
that product P, f P 2 , P3» P4 are detenlined. This leads 1 o the continuous density D ra P i c , D m P2c, 
D m rMc D m P4c for each product P u eL P 3 , P<. It is pnsible to use the already determined 
parameters p, should those already njave been detenu 1 ed in the selection step 40. 

At this point in the method the continuous densitjes D m P i c , D m i>2c, D m P3 C , D"W for each 

basis of those densities D m P i c , D'W, 
Pu P2, Pa, P* can be determined. It is 
, the found continuous density D"W, 
<ct the normalized discrete demand 



product P,, P 2 , P 3 , P4 have been determined. On th< 
D''^, D m P4c the safety stock level Si for all the produc 
however possible that due to algori|imic uncertainty s 
D nl P2c, D" 1 ,.^, E>"V4c does not accurately enough ref 



density D N p ld , D N PM , D N P4 d. A 



typical case in wjich this could happen is one where a 



has been determined for each cluster 
ongs to the respective cluster Q, C 2 
demand density DVd, D n p3j, 
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final distance dc(X>, D N Pi d ) between the cluster rk oresentative Xi, X 2 and a normalized 



discrete demand density D N n d » D^L D K wd, l^V* thfi 
to be a member of the corresponding cluster C u C 2 , 
the clustering step 30 is to perfonr 
to a random pairing selection, rrbaning selecting 
densities D N Pld7 D N p2d> D N p3 dy D^U one after tto 
randomizing algorithm and also oil the number of 
the above described effect of a large final distance 



found continuous density I^Vic, ufc*P2c> D^wc, D m F 
density D N Pld , D N wd, D N p3 d » D*W 

i 



13 
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is in the clustering step 30 determined 
relatively large. A typical scenario for 
the determinaticU of the distance dc(Xi, D N Pid ) according 
rom all normalized discrete demand 
other at random. Depending on the 
ations, situations can occur in which 
(Xj, D K rid ) may be present. This large 



final distance dc(Xi t D'Vd) would (eventually lead [id an undesired distance between the 



Sbe used to reduc 
| continuous 



A clustering correction step 60 can 
?3> P4 a distance d p between its 
normalized discrete demand density D*wd, D"™, D"] 
a threshold test step 62. For those c'f the products Pi 
predetermined threshold value tb, I a substitute-sel^ 
condition that the threshold value th is exceeded, 
D N pi d , D N !>2d, D N P3d> D N P4d is enabled, jdepicted in for 
to the substitute-selection step 63. ! 

In the substitute-selection step 63 fcjr the correspom 
D'Vid. E> N i>2d, D N W> D N P4< t out of the 1 ^determined sei 
a substitute continuous density D^iL ©"whs, I^Ws, 
that product Pi, P?, P 3 , P 4 as its ccintinuous density 
words, instead of using the cluster-representative 
selection step 40, the selection step 40 is again 
normalized discrete demand density j D*™, :D n pm, D* 
substitute continuous density I^ricsjl^j^cSvErpacs* D 



form of a logical OR-function 66 iii fig. 1, for the 

E^wc, that is only enabled, depicted I in form of a logical ANEMuncuon 65 in fig. 1, to be 
used if the threshold test step 62 resulted in the fin ling that the threshold value t* is not 



exceeded. 
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and the normalized discrete demand 

this problem. For each product Pi, P 2 „ 
D"Vic, EPpzc, D m P3c, EPimc and its 
I, I^Vw, is determined. This is done in 
?i> P?, P* whose distance d p exceeds a 
lion step 63 is performed. On the 
normalized discrete demand density 
< if a logical AND-function 64, as input 

normalized discrete demand density 
of continuous model densities D°mc 
Ws is selected, which is defined for 
5 m P,c, D m P2c» D m P3c> D"W. With other 
screte density R rc as input to the 
jerforjned using now as input the 
d. D N ?4d itself. The found normalized 
4cs is used as a substitute, depicted in 
: mtinuous density D m PIC> 
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Once the continuous densities D m r&, D^ci D^cDW have been determined, an adjustment 



step 70 follows, wherein for each product P,, P 2 , P * 



D m rec, D m P4c is adjusted into the conjtimious demand ;iensity Dpi C) Dp2c, Dp3c, Dp*, such that the 



continuous demand density Dpu, pnc> Djw c > Dn* 
Dp2d, I>P3d, Djm, respectively, havej substantially ths 
Therefor the factors s 0 , s v from tfce norriaalizatior 
adjustment factor set v 6Si , v^, v^, v,^ iare used 
rescaling. 



Thereafter for each product Pi, P 2 , P 3) P 4 the corresp rtiding continuous demand density D PJC 



10 Dpao Dp3 C , D P4c has been determined. These continubus demand densities Dpi c , Dp2c, Dp3» 



5 
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P 4 the continuous density D m P io D m P2 C , 



;tiid the discrete demand densities D n< i, 
same average and the same variance, 

5tep 20 are used, more precisely, the 

i 

jThe effect is hence a reshifting and 



with 

X 



with the lead time LT for the specific product Ps 
demand density Dpic of that product Pi, the probab 
product Pi being demanded during t|e lead jtime LT 
of k days, or more generally time u^its used in the cfetnand time series H(Pi), H(P 2 ), H(P3>, 
H(P4>, more than 0 pieces of that product Pjlare being demanded during the lead time LT, 



safety stock level s* for each of the 



Dp4e can be used for determining 1 ) a corresponds, 
products P, f P 2 , Pa, P* in a safety] stock detenhimtion step. In addition herefor, a given 
service level SL, that here is an ai-ijervice level, is u ;?d. 

A possible formula for determining the safety stock 



elvel Si is the following: 



and the average p. of the continuous 
lity po of demand of 0 pieces of that 
a ad the probability pk that on a number 



Generally speaking, the continuous Idemand density determines for the specific product 

5 ill 

Pj the relation between the given service level SL jof t ijat product Pi and its safety stock level 

Si. Having in addition for that specific product Pijthe 

f ) ; : 

the time that passes between placing an order foir fcj : 



and the factual delivery of the product P t toftie onlerlnjg entity, the minimum quantity Q rn of 
the product Pi in order to be able to luarantee thejgiv sii service level SL can be determined. 



corresponding lead time LT, which is 
that product Pi to a product provider 
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determi ed 



si ck 



•i;- is 

: 

The minimum quantity Q m (P») is then given by 
product Pj during the lead time LTf and 

Pi. This also means that as soon is the currten;t 

v 

below or reaches that minimum quantity <J m , a 
the reorder point. The reordering jvill than goarantfefc 
S(Pi) will not be too low, A reorcj aing quantity Q 
minimum quantity Q m such that tiie reorder prdcei 
can be held low. Since the reordering quantity 



reordering costs, warehousing costL depreciation, 

i t . i 
stock controller to determine the 4 ppropriate 



reor aing 



Qr( ) 



e 



Julc 



i; 



irebrd ing 



reordering, therefore the system wduld be ^provided 

reordering quantity Qr(P«). However it is also pdssil^l^ 

it : * 

predetermined to enable an automatic reo: 



sum of the average demand for that 
safety stock level S| for that product 
level S(P,) for this product P f falls 
should be initiated This is called 
that the product's current stock level 
?0 should be typically higher than the 
and possible associated reorder costs 
depends on associated costs such as 
it can be put at the discretion of the 
quantity Qr(Pi)- For the automatic 
y the stock controller with the desired 
to have the reordering quantity Q r (P|) 
without external intervention- 



pre c ess 



There are various possible definitions of service 'itevi 
defined as the percentage of cycle ({times Without 
i.e. the inability to satisfy an occulting derharid for 
policy can be applied that triggers pr initiates^ 



SL, and herein the service level SL is 
ajW occurrence of a stock-out condition, 
product Pi on stock. A replenishment 
jproppct order as soon as the current stock 
level S(Pi) thereof drops to or belodf the mi&nium qx t ntity Q rt 



i 



safety 
wice an order 




Hence, once having determined 
performed wherein via an output 
towards a product-providing entity! in the:ev£nt thftl 
product P,, P 2 , P 5 , P 4 falls below 
automated way of controlling a 
level Si- 



The precision of the above desc ibed riietfaoe 
exceptional rules. The value of 0 j ieces cjf pioi 



miniinum quantity 



ict stock level 



a 11 



be improved by applying further 
Pi, P 2 , P3, P4 demanded is insofar a 
special value since it is in some s tnse an in^gpifi rant value. This anomaly of selective 
significance can be taken into accou it by tb# fi$l<£wij ^ improvements: 

The demand for 0 pieces of the products Pi/j^t, 



In the normalization step 20 the 



deiband va&uejf 



the discrete demand densities D PId , 



IW pfU i* 
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UC!> 



^vel Si, a product order step can be 
a product P lr P2, P 3 , P 4 is initiated 
the current stock level S(P0 of that 
Q m . Thus, the method provides an 
rfider use of a determined safety stock 



P 4 i i assigned a separate bin. 

forO fctiec s of the products Pk P2, Pj, P4 in 
! nversely reseated with respect to the 
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demand values for more than 0 pijbces of ^he same product P Jt P2, P3, P*> This can be done 
by detennining the average zero-i emandijprobability p a o for the demand for 0 pieces of a 
product Pi. P 2 , Ps, Pa, over all produces -P Jt P 2 , P 3 , P 4 , and determining the average 
non-zero-demand probability pu+ fpr the dbmand for more than 0 pieces of a product P l7 P 2 , 
P 3 , P«, over all products P,, P 2 , P 3 , f * Furt|ermore a : elation level l r is defined that gives the 
desired relation between the average nof-zero~den and probability p n i+ and the average 
zero-demand probability p a0 * Then the resting is executed in that the demand values for 0 
pieces of the products Pu P*t Bk are e^ph multiplied by a zero-demand rescaling factor 
which is the relation level l r dividpd by &e average 
demand values for more than 0 piqpes of t|e product! 
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non-zero-demand rescaling factor which is :1 minus 
average non-zero-demand probabili y p a i+, J 1 



This rescaling has the effect that flbr slowimbving ii sms the true demand occurrences are 
not outscaled and therefore reduce| in their significa ice by the demand value for 0 pieces. 
Otherwise all slow-moving items vjfould re|ult to am 
density for more than 0 pieces nfay difffr istronglj 



significant values, i.e. those above ft pieces 
significant, it should not be suppressed hi 
result of the clustering. In the casejfof sup 



zero-demand probability p a0 , and the 
Pu P2. P3. P* are each multiplied by a 
the relation level l r divided by the 



>e in the same cluster Q> although the 
The rescaling allows to stress the 
the demand value for 0 pieces is less 
A suppression would again distort the 



Although 
since 

re^sion, sfow-moving items would be clustered 
together with non-slow-moving iteifis, whiJh would a so not be correct, and even within the 
class of slow-moving items there jfe a distance between "slow" and "very slow", i.e. the 
products Pi, P2. Ps, P4 with very few demaid occurra ices, e.g. 3 times in a year, and those 

with more, e.g. 20 times in a year. jThese dtstknees shall not be neglected in the clustering 

1 ill 

step 30, hence the rescaling is seen is the *>propnate 
in a more accurate clustering in part cular frftr the slow 



In the selection step 40 the 



measure. The rescaling hence results 
moving items. 



value for 0 pie ces of the products P x> P 2 , P 5 , P 4 is 



suppressed. Again this is particidafly useml because 



model densities D^mc the value of CB pieces! 



model density D m M e. Now here a |ptal suppression 
30 clustering has already been done. 



0uld disi on the assignment of the best fitting 
las no detrimental effect, since the 
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In the adjustment step 70 the deim nd vahi j for 0 pieces of the products P,, P 2 , P*, P* can be 
reused and iescaled. This achieves an uncj storted eq rivalence between the original discrete 
demand densities Dm, Dp2 dJ Drm, ,! >i*w an<! the final continuous demand densities Dp, c , Dp2c, 

The described method can be exe nxted c i 4 compt iter that is equipped with a computer 
program that when executed car tes oui the steps of the method. More precisely the 
computer would be designed to © inprise a process >r for executing the different steps, a 
memory for at least intermedial ly stor ngj the re suits of the different steps and the 

i l 

predetermined set Ld of continu rus mc ie| densii ies D*W The computer could also 
comprise an input device for recer ing frcj a it user a c 

and/or the reordering quantity Q* The c aniputer c ould furthermore be connected to a 
computer network and electronical y effec tile issua ice of an order to a product provider 
for a specific product Pi, P 2 , P 3 , P* it a poi it in time *hen the current stock level Si of that 




r 6r. In th|s sense the method would provide an 



product P,, P 2 , P 3 * P4 falls below th i 
a predetermined reordering quantir Q m he 

automated solution for the problen of maintaining rife stock of several products P,, P 2 , P 3 , 
P 4 at a level that corresponds to a di sired sc vibe level SL. The automatic order issuance can 
also be effected by controlling an o ^tput dj vice of an j kind, e.g. printer, plotter, telephone, 

etc- wl osp outpu 



fax machine, loudspeaker, display, 
e.g. issuing an automatic order pric but th* 



an order e-mail that is automatic illy sei t to the 



over a phone connection a voice me isage t^ it is receh able at that product provider, creating 



computer could also instead of or ering i it&ctiy su< :h a product P u P2* Pa, P* issue an 



alerting message to a stock control! r who i tiepeby haj 



decide thereupon whether or how 
ordered. 



nuch 0 



In fig. 2 a product stock control system i 
method. The system comprises a 
P 3 , P4 its demand rime series Hi 
LT(P^), LT(P 4 ), its reorder quantity 
S(P,), S(P 2 ), S(P 3 ), S(P 4 ). It 
product-specific service level S 
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infd mint num quantity Q„,. The system can use 



is receivable by a product provider, 
can be s|bt to the product provider, or issuing 



iroduct provider, and the like. The 



the resp< etive product P u P 2 , P3, P* is to be 



djepicted 



H(pJ, |H(P3), 
Pi), 



SL(P 3 ) 



a stock level control system and can 



making use of the above described 
uterp tjiiat recAves as input for each product P,, p 2 , 
l(P,)> its lead time LT(P,), LT(P 2 ), 
VCf 2). QrC S)> Qr(Pd) and the current stock level 
receive! a service level, that can be a 
SL(Pd). The computer executes by 
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means of its processor the stoi *d con pujter pro pram that performs the safety stock 



determination method, determii ing g< 
normalization step 20, clustering s ep 30, 
denormalization step 70 for each j noduct 



Dp2c, Dpac Dp4c is determined- Thei \ follow >.the safety 



for each product P,, P 2 , P 3 , P* the 



in dependence on the product-spec: Be serv ce levels: .(Pi), SL(P 2 )» SLCPj), SL(P 4 ) 



In a subsequent reorder-point det irminaiflooj step 9|3 ; 
Qm(P 2 ). Qr»(p3>, Qm(P<) is determia *d for 
10 times LT(P,), LT(P 2 ), LT(P*), LT(P J). 

In a subsequent stock analysis stej ! 100 ufid^r ^se c f 



P 3 , P 4 . As soon as the stock situatic n for al leiast one 



mg e.g. 

i 

election s 
Jp*. Pa, 
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lafety s ocjc level 



through the conversion step 10, 
£p 40, parameter selection step 50 and 
the continuous demand density D P i c , 
stock determination step 80 in which 
fc(Pi), Si(P 2 ), $i(P 3 ), Si(P 4 ) is determined 



i.e. the minimum quantity Q m (Pj), 
h produbt Pi, P^ Pa* Pa under use of the lead 



the current stock level S(P,), S(P 2 ), 



S(P 3 ), S(P<) the computer 1 detent ines th| current s ock situation for each product P u P 2> 



Of the products Pi, P 2 , P3, Pa becomes 



critical, i.e. the current stock lev* 1 S(Pil S(P£), SpP 3 ) 7 S(P 4 ) approaches the minimum 
quantity Q m (Pi), Qm(Pj)t QmCPa), Q kP*) tcj a predetermined extent, e.g. becomes equal or 
lower than the minimum quantity ^Q m (Pil : QmCP2)» 
perform a stock control action for hat prc|iu0t £w Pj, 



possibility for the stock control acti in is tofcerfoirn a 
automatic order 102 under use of t te conwpcpnding 

Qr(P4). 



It is obvious for the person skilled in the lit |tfa& the 



3m(Ps), Q*(P*), the computer 1 will 
P 3l P4. The stock control action can 



e.g. be a stock-alert signal 101 t kat is lent put vGa a computer output port. Another 



• i 



hardware, software, or a combinatic 1 of those/ Also, 

fashion on one single computer sysi sm, or Li d&tribi 

' |. 1 ; 1 :[ 

25 are spread across several interconnected coppptep or 



of a computer system - or other ap; aratus adapted fop carrying out the methods described 
herein - is suited. A typical comb hationfcf 
purpose computer system with a co: hputerpro; 



gram 

controls the computer system suchtthat it tames ou 
present invention can also be embe ded inla comput* 
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jroduct order step 1 10 by initiating an 
i border quantity Qr(Pi), Q,(P 2 ), Qr(P:0, 



.present invention can be realized in 
can be implemented in a centralized 
ted fashion where different elements 
computer systems, whereby any kind 



hardwgre and software could be a general 



ttjfat, when being loaded and executed, 
. the methods described herein. The 
€ f program product, which comprises 
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all the features enabling the impli mentat; >n of the 
when loaded in a computer system j- is abl 



Computer program means or comp ater pre gram in thi 

■ '. i 

i ji 

seti 6f instructions 



in any language, code or notation, M a 
an information processing capabiii iy to pe form a 
either or both of the following a | conve 'iion to 
reproduction in a different material jfomu 



Any disclosed embodiment may b 



shown and/or described. This is als ipossil le for one 



-19- 



pethods described herein, and which 
to cany cut these methods. 



particular 
another 



>mbii ed with one 



10 It is obvious that a person skilled: la the can 



ways without departing from the 
subsequent claims. 
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fe or several of the other embodiments 
?r more features of the embodiments. 



mollify the shown arrangements in many 
the invention which is encompassed by the 



! ! 
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present context mean any expression, 
intended to cause a system having 
function either directly or after 
language, code or notation; b) 
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i 1 mi ' 

1. Method for assigning to a set of products (P u P2, Pj, P4) a set of corresponding 

il . 1 J l» 
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continuous demand densities fpnl Dj, if Dp*., Drac) comprising the following steps: 



(H(Pi), H(P 2 ), H(P 3 ), H(P 

Dp3 d> Dp4d) 7 




- a conversion step (10) wherein forleaclijprdduct (P u P2, P*, P4) its demand time series 

cbnv^i fed into a discrete demand density (Dpu, Dpza, 



said 



- a normalization step (20) whe 

. H . it 

are transformed into normalized discr^l e 



1 1 

< iscrste demand densities (Dpi*, Dp**, Dpsd, Dp4d) 
demand densities (D N P id, D'W D n p*u D N P4d), 



- a clustering step (30) wherein e&ch of sja 
D N p2d» E^pw, D N wd) is assigned to a cjt 
(Ci, Ca) a cluster-representativl discretjs 



a selection step (40) wherein f< 
a predetermined set (Ld) of o 
continuous density (D" cc) is 




: normalized discrete demand densities (I^W 
ter (Ci. C 2 ) and wherein for each said cluster 



lenity (Rfo, Rra) is determined. 



representative discrete density (Rjc) out of 
)del densities (D m Mc ) a cluster-representative 



- a parameter-determination sted (50) wj oeig for each product (Pw P 2 > P3, P*) for its 
cluster-representative continudps dpnsii KP^cc) product-individual density parameters 
(p(P)) are determined under uisfe of whii fauFbr, each product (Pi, P 2 , P?, P4) a continuous 



density (D^ic, D'W, D m w c , DiL) i$ deti 



a adjustment step (70) wherein' lor each 
said continuous density (D n Vi c ;|D m i>ic, 
demand density (D P i c , D^* D?4 Dric). 



ed, 



ri^dact (P w P 2 , P 3 , P«) average and variance of 

10 ! 

^ D m Pte) are adjusted to form said continuous 



2. Method according to claim 1 .ftdther 
wherein for each product (Pi, 1 1, P3, P4; 
(D m pic» D m ?2o, D™*, D ra P4 c) and its hoi 
D N P3d, D w P4d), is determined aj^d whe 
whose distance (dp) exceeds 



>n prising a clustering correction step (60) 
[distance (d p ) between its continuous density 
iz^d discrete demand density (D N P id, D N p 2 <i, 
J ^fiir those of said products (Pi, Pi, P 3 > P4) 
j&rmined threshold value (t h ), in a 
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10 



15 



20 



substitute-selection step (63 > for t 
density (D*™, D N ™, D*™, 1)%) o 
model densities (E^mc) a substitute co 



selected, which is defined f< 
(EF Plc , D^, D m P3 C , D m P4c). 



that p 




3. Method according to claim t or 2, ^ 
time series is substantially se| negated 



4. Method according to one of claims 
normalization information is r|aintiam< 

5. Method according to one of cjfctims I f 
initial cluster representative C&u Xa) 
comprises for each of the no^nalized 



- a distance determination step whejreinj 
X2) its distance (dc(Xi, E^fw)) 
determined* and 

- a modification step wherein lie onej 
towards the normalized discre se dfci 




distance (dc(Xu D*™)), said 
defined as new initial cluster r< 



modifi 



i 

sent 



6. Method according to claim Wher 
modification step are repeated one or 
demand densities (D*™, D* ?m, ID*! 
representatives (X,, X 2 ) that arJpreientl 
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|cpiresponding normalized discrete demand 
said predetermined set (Lj>) of continuous 
aiaous density (D% lcS , D'^, D"^, E^ P4cS ) j s 
act (Pj, P 2 , P 3> p 4 ) as its continuous density 



in the conversion step (10) the demand 
^demand trend and periodic demand. 



St, wherein in the normalization step (20) 
[Of the adjustment step (70). 



^wherein for each of the clusters (Cu C 2 ) an 
iyided, and wherein the clustering step (30) 
ate demand densities (D"™, D"™, D^*, 



|t pach said initial cluster representative (X w 

i 

salized discrete demand density (£>%) is 



III 



cluster representative (XO is modified 
ity (E^Pid) to which it has the smallest 
t$ter representative (Xi) being thereafter 



3ke distance determination step and the 
times for each of the normalized discrete 
W, starting with the initial cluster 
the previous last modification step. 



1 
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7. Method according to claim 5 
that are present after the last 



or 6' 
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a- r 



, wieieini the initial cluster representatives (Xi, X 2 ) 



nodificmon step are selected in the clustering step (30) 
33 the cluster-representative discrejte dpisityj(Rfc). 

8. Method according to one of dlaiiES 1 lb 7 ,j wherein in the clustering step (30) a said 

subslkntiaily not modified when its distance (dc(Xj, 



cluster representative (X u X2) is 

j| 

E^Rd)) to those normalized dfscreie 

J 



ian(J densities (D N Pld , D N P2d, D N P3l i, D N ^) that 



belong to the corresponding cjbstei r (CKj), is below a predetermined value. 

9. Method according to one of jplaiins Jjto 8| , wherein the predetermined set (L D ) of 
continuous model densities p"\jb) Jptaihs the model densities (D w Mc ) with its 
parameters (p) in an undetennfinedlfoiJ and jwherein the parameters (p) for the clusters 
(Ci, C 2 ) are determined in the |elec tionljtep j40). 

10. Method according to one of aainjs 1 fi 9 > Wherein in the normalization step (20) the 

[ I ! 

demand value for 0 pieces 01 the prlpucts (Pu P2, P3, P4) in the discrete demand 
densities (Dpu, Dp**, is resfiled inversely to the demand values for more 

than 0 pieces of the same prodlct P 3 , jP 4 ). 

11. Method according to one of|clains || to IjO, wherein in the selection step (40) the 
demand value for 0 pieces of die pi aduKs (Pi P 2 , JP 3 , P 4 ) is suppressed 



12. Method according to one of cl 



1 to 11, further comprising a safety stock 



detenninaiion step (80) wherean v ithlgiv^i service level (SL(P,), SL(P 2 ), SL(P 3 ), 
SL(P 4 )) from the continuous dejmar d d&sityj(Dpi C , Dpto D P3c , EW) for one or more of 
the products (Pj, P 2 , P 3 , P4) a prodiict ftfetylistodc level (s,(P0, Si(P 2 ), Si(Pj), Si(P«)) is 



determined. 



13, Method according to claim 12 



fuimerBoinp rising a product order step (110) wherein 



via an output device an order 1 or a 
product-providing entity, in thj 
S(P 4 ) of that product (P u P 2 , P 



, P^ 

Qm(P3)> Qm(P4)) that is the sun of the Ijvera je of the discrete demand density (D^ Id7 
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sai* product (P,, P 2 , P 3 , P 4 ) is initiated towards a 
e\Jent Biat c| current stock level S(P|), S(P 2 \ S(P 3 ), 
falR bel jw a minimum quantity (Qm(Pi), QmCP*)* 
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Dp**, Dwd, Dpad) for said prodtat (P u A P 3 , P 4 ) <Wing the product's lead time (LT(P,), 



LT(P 2 ), LT(P 3 ), LT(P,)), and 
Si(P<0) for said product (Pi, P 2 , 



15. Computer program product 



means stored on a computer-readable medium. 



_ ^ j 
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the detjmiaed safety stock level (Si(Pi), Si(P 2 ), Si(P 3 ), 
Pi, P«). 



14. Computer program product q Dmprisiig program code means for performing a method 
according to one of the preceding cla 



accc 
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to claim 14, comprising the program code 
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scr 
sc 



The invention is directed to a method 
corresponding continuous demand deisi 
demand time series is converted into a fis 
discrete demand densities are transforn^ed 
clustering step each of the normalized dp: 
for each cluster a cluster-representative Si 
each cluster-representative discrete den|ity I 
densities a cluster-representative { cd 

| i. 

parameter-determination step for F ac ^ P 1 
density product-individual density paratoet^ 

product a continuous density is dlten^inej 

El i! 

continuous density is adjusted into tjpe c^pnti 

j. . 

\ 



tACT 

ii 

)r assigning to a set of products a set of 
L In a conversion step for each product its 
ke ciemand density. In a normalization step the 
ato normalized discrete demand densities. In a 

I is 

I P 

|e demand densities is assigned to a cluster and 

v. 

te density is determined. In a selection step for 
it cjif a predetermined set of continuous model 

ItiMous density is selected. In a 

I Ij 

[tact! for its cluster-representative continuous 

ii 

aiie determined under use of which for each 
Inji an adjustment step for each product the 
jiouS demand density. 
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